Stripes and holes in a two-dimensional model of spinless fermions and hardcore bosons 
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We consider a Hubbard-like model of strongly-interacting spinless fermions and hardcore bosons 
on a square lattice, such that nearest neighbor occupation is forbidden. Stripes (lines of holes across 
the lattice forming antiphase walls between ordered domains) are a favorable way to dope this 
system below half-filling. The problem of a single stripe can be mapped to a spin- 1/2 chain, which 
allows understanding of its elementary excitations and calculation of the stripe's effective mass for 
transverse vibrations. Using Lanczos exact diagonalization, we investigate the excitation gap and 
dispersion of a hole on a stripe, and the interaction of two holes. We also study the interaction of 
two, three, and four stripes, finding that they repel, and the interaction energy decays with stripe 
separation as if they are hardcore particles moving in one (transverse) direction. To determine the 
stability of an array of stripes against phase separation into particle-rich phase and hole-rich liquid, 
we evaluate the liquid's equation of state, finding the stripe-array is not stable for bosons but is 
possibly stable for fermions. 
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I. INTRODUCTION 



Stripes have been an area of active study in the high- 
temperature superconductivity research. They are mod- 
ulations of charge and spin densities, and can be static or 
dynamic. In 1995, Tranquada and coworkeraj observed 
the coexistence of superconducting and stripe domain 
order for the material Lai.6-2;Ndo,4Sra;Cu04 at around 
doping X = 0.125. The stripes are one-dimensional ob- 
jects on a two-dimensional plane, and thew have been 
called "self-organized one dimensionality. "□ Search for 
stripes in the Hubbard and t — J models has been intense 
and is still going strong (see Ref. ^ for a review of some 
experimental results and theoretical considerations). 



^From the theory side, there are basically two roads. 
One approach is using microscopic models such as the 
Hubbard and t — J models to study numerically stripe 
formation and phase separation. DMRG appears to be 
the best tool in such studies and has been applied to the 
t — J model in Refs. ^,|[ Another approach to stripes is 
to use effective field theories, treating stripes as macro- 
scopic objects (macrosopic in the sense of many-particle 
collective motion). For example in Ref. |^, Zaanen con- 
sidered stripes as "a gas of elastic quantum strings in 
2 + 1 dimensions." The microscopic, numerical approach 
can be used to understand the mechanism for stripe for- 
mation and to derive macroscopic parameters such as 
stripe-stripe interaction energy and stripe stiffness, but 
it is limited by computational power. The macroscopic 
approach on the other hand takes these macroscopic pa- 
rameters as an input and can obtain results relevant to 
experiments, but it can explain less about mechanism. 
We will follow the microscopic route which enables us 
to obtain macroscopic parameters from diagonalization 
data. 



A. The model 

In this paper, we study stripes in a strongly-interacting 
model of spinless fermions and hardcore bosons on the 
square latttice. The Hamiltonian is 



H = -t 



£( 



V 



(1.1) 



with periodic boundary conditions. We study both the 
spinless fermion and the hardcore boson versions of the 
model. c| and Ci are the spinless fermion or hardcore bo- 
son creation and annihilation operators at site i, fii — clci 
the number operator, t the nearest-neighbor hopping am- 
plitude, and V the nearest-neighbor interaction. At the 
each site there can be at the most one particle. Further- 
more, we study the strong-correlation limit of the model 
and take V = +oo, i.e., infinite nearest-neighbor repul- 
sion. 

The spinless fermion and hardcore boson model with 
infinite nearest-neighbor repulsion involves a significant 
reduction of the Hilbert space as compared to the Hub- 
bard model. The Hubbard model on the 4x4 lattice 
at half-filling, with 8 spin-up and 8 spin-down electrons 
has 1,310,242 states in the largest matrix block, after 
reduction by particle conservation, translation, and the 
symmetries of the square J3 In our model with infinite V, 
after using particle conservation and translation symme- 
try (but not point group symmetry), the largest matrix 
for the 7x7 system has 1, 906, 532 states (for 11 parti- 
cles). We can therefore compute for all fillings the 7x7 
system whereas for the Hubbard model 4 x 4 is basically 
the limit. This also means that at certain limits we can 
obtain results that are difhcult to obtain in the Hubbard 
model, for example, we can exactly diagonalize a system 
of two lcngth-8 stripes on the 8 x 28 lattice. 

This is one of the two papers that we are publishing 
to systematically study the phase diagram of the spinless 
fermion and hardcore boson model on the square lattice 
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with-infinite nearest-neighbor repulsion. In the othe r pa - 
per,Q we study the dilute limit of this model Eq. (1.1), 
focusing on the problem of a few particles, and that limit 
is dominated by two-body interactions. In the present 
paper, we study the near-half-filled limit of our model, 
where stripes are natural objects. 

With infinite nearest-neighbor repulsion, the maxi- 
mum filling fraction (particle per site) for our model is 
1/2, for both spinless fermions and hardcore bosons. At 
exactly half-filling, the particles form a checkerboard con- 
figuration that cannot move. Adding a single hole to it 
cannot produce motion either because the hole is confined 
by neighboring particles. A natural way to add holes to 
this system is to align them in a row going across the 
system, as indicated in Fig. |^. We call this row of holes 
a stripe and it is the subject of study in this paper. 
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manganitesjii the formation of inhomogeneities is domi- 
nated by lattice distortions and orbital degeneracies (not 
to mention static pinning disorder), so that quantum- 
fluctuating defects find no role in current theories of man- 
ganites. (The-aame may be true for the Verwey transition 
in magnetiteO which was modeled previously with spin- 
less electronstj. ) Conceivably the spinless model could 
be realized in an apjay of semiconductor quantum dots 
in a magnetic field. ta 

The most plausible realization of the boson version of 
our model would be an adsorbed gas of ^He on a sub- 
strate with square symmetry. (The same model - but on 
a triangular lattice - has beecuintroduced independently 
to model ''He on nanotubeEjcd) . The fermion case could 
similarly be realized, in principle, by spin-polarized "^He, 
however it would be difficult to achieve a degree of po- 
larization sufficient to neglect the minority spin state. 

This paper appears to be the first using exact diagonal- 
ization (as opposed to DMRG) to investigate the prop- 
erties-of interacting stripes with a microscopic Hamilto- 
nian.y 



B. Stripes in spinfull models 



FIG. 1; An isolated hole (shown in the box) in a half-filled 
region cannot move. A row of holes (a stripe) can slide along 
in the vertical direction, with the arrows showing the possible 
moves. 

For many systems with certain aspect ratios (e.g., the 
6x7 lattice depicted in Fig. |l|), the stripe state is the one 
with the smallest number of holes and nonzero energy. 
Here we have an interesting kind of physics of extended, 
fluctuating, quantum mechanical objects, that are results 
of collective motions of many particles. 

Spinless models related to the present one have been 
invoked occasionally in the literaturdj but usually in the 
context of a specific question about spinfull models: it 
was recognized that a spinless model would capture the 
same physics with fewer complications. The only sys- 
tematic-jwork on phase diagram of spinless fermions is by 
UhrigJ13 However, the method is expansion around infi- 
nite dimensionality; this will have special difficulty with 
arrays of domain walls, since wall fluctuations and the 
attendant kinetic energy are strongly suppressed in high 
dimensions. 

The Falicov-Kimball model is an alternate way to sim- 
plify the Hubbard modcl,_which is also commonly pre- 
sented in a spinless forBi,iij and develops stripe- like in- 
commensurate patterns.li3 However, this model includes 
a second immobile (classical) species of electron, so its 
stripes cannot have quantum fluctuationsJij 

No known electron system realizes our madaLeven ap- 
proximately. A "half-metallic" ferromagnetcfltJ - mean- 
ing that for one spin state, the conduction band is all 
filled or all empty - realizes a spinless model for the other 
spin state. However in the best-known half-metals, the 



In real fermion systems, stripes (antiphase domain 
walls containing holes) have been a promincaL. olpj ftct 
of experimental study in the high- Tc cuprates.BQcSlHj-Eil 
Stripe-based mechanisms have|-|bLe|en proposed for high- 
temperature superconductivityJ3'Ej but the prevalent cur- 
rent gpi nion, is that stripes compete with supcrconduc- 
tivityH'BEJ 

Nevertheless, stripes are obviously clues to how the 
charge and spin degrees of freedom interact with each 
other, which is important in a majority of the high-Tc 
theories. Stripes are modeled theoretically and numeri- 
cally using the same Hubbard model (or its variations) 
which were already accepted as models of homogeneous 
phases. It is still unsettled whether stripes are stabilized 
in this model, which omits long-range and even nearest- 
neighbor Coulomb repulsion. The most explicit calcu- 
lations are by Density-Matrix Renormalization Group 
(DMRG), .adapted to two-dimensional systems formed 
into strips D'EI The results favor stripes with 1/2 holes per 
unit length, as found in experiments; however, those sim- 
ulations are strongly influenced by the (necessary) open 
boundary conditions on two sides. 

What relation can our spinless model have to these 
spinfull models? Firstly, at the intermediate scales ap- 
propriate to a stripe array, the spinfull and spinless sys- 
tems are modeled in very similar ways: stripes are mu- 
tually repelling, quantum-fluctuating strings. This level 
is appropriate to studies of the anisotropic transport ex- 
pected in a (non-superconducting) stripe array, as well 
as the,stjiac interactions, as has been emphasized by Za- 
anen.l3'E3'EZl In the spinfull models, only a little has been 
done to explicitly relate the macrpscopic parameters to 
the underlying microscopic model E3 In our model - and 
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unlike any spinfuU model - exact diagonalization can ad- 
dress phenomena such as the effective mass and interac- 
tions of carriers on a stripe, long- wavelength stripe fluctu- 
ations, stripe-stripe repulsion, or the transfer of carriers 
from one stripe to the next. Thus, not only analytically 
but numerically, one can go quite far in extracting such 
parameters for our model. 

Secondly, if spinfuU stripes are stabilized at all (in the 
absence of long-range repulsion), it is by kinetic energy 
in the same fashion that our stripes are stabilized. In ei- 
ther fermion belonging to one domain hops trans- 
verse to the stripe, and (since the stripe is also an an- 
tiphase domain wall) still finds itself correctly placed in 
the order of the new domain. In our model, that order 
is the "checkerboard" pattern, which is enforced by the 
nearest-neighbor repulsion V. From this viewpoint, V is 
analogous to the antiferromagnetic coupling J in the t-J 
model. (Our V oo limit then corresponds to the case 
J ^ t, which is the unphysical regime of the t-J model.) 

An even more precise-iprrespondence would be to 
stripes in the t-Jz model:caE3 in that system (as in ours) 
the undoped, ordered state breaks a discrete (Ising-type) 
symmetry, and hence is practically inert, supporting no 
gapless Goldstone (spinwave) excitations: only the stripe 
itself has low-energy excitations. If the ordered state had 
a continuous symmetry, as in the Hubbard or t-Jjsaod- 
els, spin- wave modes can mediate a, attractiorH be- 
tween stripes (where d is the stripe separation) . If short- 
distance kinetic energy were to favor a stripe array, the 
long-range force mediated by continuous spins inaplies 
phase separation in the limit of very small doping.E3 

Our stripes with an occupation of 0.5 hole per unit 
length are insulating and correspond to insulating stripes 
of 1 hole per unit length in the t-J model (the factor of 
two reflects the density of the cpieted background state), 
as implied in the original ideasESEjO about stripe stabi- 
lization from a strong-coupling viewpoint. 



C. Paper organization 

This paper is organized as follows. First (in Sec. ||) we 
describe briefly our exact diagonalization code for study- 
ing the near-half-filled limit of our model. We use trans- 
lation symmetry to block diagonalize the Hamiltonian 
matrix. A graph viewpoint motivates a method of build- 
ing the basis states from a starting configuration and 
explains the relation between boson and fermion energy 
spectra. 



Next (Sec. HI) we study the problem of a single stripe. 



Here particles can only move in one direction, effectively 
reducing the 2D problem to a ID one. This one-stripe 
problem maps exactly to a spin-1/2 chain, which is ex- 
actly solvable using Bethe ansatz. Using the single-stripe 
energy dispersion relation along the direction perpendic- 
ular to the stripe, we obtain stripe effective mass for mo- 
tion in that direction. 

In the next two sections, we study the problem of one 



or two holes on a stripe. (For the one-hole case, the 
bosons and fermions have the same energy spectrum.) 
We determine the dispersion relationship (energy gap and 
effective mass) for one hole moving on a stripe, and study 
the binding of two holes. It is informative to study the 
finite-size dependence of the energy on the lattice dimen- 
sions in both directions: in particular, it decays expo- 
nentially as a function of the lattice size perpendicular 
to the stripe, which we explain quantitatively in terms 
of the stripe tunneling through a barrier between two 
potential wells. 

One of the motivations for studying our model is our 
interest in a simple model of interacting stripes. In 
Sec. VI we have exactly diagonalized systems with two, 
three, and four stripes, and find that the stripes repel. 
The interaction energies scale as the inverse square of 
the stripe spacing, which (like the above-mentioned tun- 
neling) follows from the one-dimensional nature of the 
stripe's transver se m otion in our systems. 

Finally (Sec. VH) we discuss the stability of an ar- 
ray of stripes by fitting the diagonalization results in the 
intermediate filling limit and using a Maxwell construc- 
tion. Our interest is whether at the intermediate-filling 
limit we have the stripe-array case or a phase-separated 
case with hole-rich regions and particle rich-regions. The 
conclusion is, interestingly, that the boson stripe- array is 
not stable and the fermion stripe-array is very close to 
the stability limit and is possibly stable. 

Our earlier publication, Ref. |3^, contains some of the 
results of this paper in a condensed form, with a focus 
on the stability of the stripe-array. The present paper 
contains more systematic and more updated results and 
a substantial number of important new results on, for ex- 
ample, boson and fermion statistics, stripe effective mass 
and stripe-stripe interaction. 



II. BASIS STATES AND DIAGONALIZATION 
PROGRAM 

In this section, we describe our exact diagonalization 
code. We begin by describing the conventions for labeling 
basis states, which are crucial to keep track of sign factors 
in the fermion case. Then we introduce a geometrical way 
to view the basis states as nodes on a graph, and apply 
this to classify the conditions under which the boson and 
fermion problems have the same. 



A. Basis states 

Each basis state with M particles corresponds 1-to- 
1 to a configuration, which is array of the M occupied 
site numbers (ji,i2, ...,H/)- Any configuration with two 
nearest neighbors occupied is excluded because of the 
infinite V in the Hamiltonian Eq. (1.1). Periodic bound- 



ary conditions are specified by two lattice vectors Ri 
and R2 (which will be called boundary vectors in the 
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FIG. 2: Square lattices with periodic boundaries: (4, 0) x 
(0, 5) on the left and (4, 1) x (1, 5) on the right. Site numbers 
are shown, following the numbering convention, upward and 
right ward. 



such that for any lattice vector r we have 
f rt2R2 = r, where ni and rt2 are two inte- 



following) , 
r + niRi 
gers. 

The configuration by itself is insufficient to define the 
basis ket, since we need to specify its sign or phase factor. 
To fix the fermion sign, we must estabUsh an arbitrary 
ordering of sites and always configurations in this order. 



Then we denote the basis state 



^11 12- 



|0>, 



The site ordering convention used in this paper is to 
start with site at the lower left corner and move upward 
along a lattice column until we encounter the cell bound- 
aries defined by the boundary vectors Ri and R2, then 
shift rightward and repeating in the next column, pro- 
gressing until all sites have been numbered. Fig. shows 
two example systems: (4, 0) x (0, 5) and (4, 1) x (1, 5). 

In this basis the Hamiltonian Eq. (1.1) acts as 



H\n) 



-t) 



s(n,m)\m), 



(2.1) 



where Ai denotes the set of states created by hopping one 
particle in \n) to an allowed nearest-neighbor site. For 
bosons s(n,m) = 1 always and for fermions s{n,m) — 
±1. And the matrix element is (r7i|_ff|n) = —s{n,m)t if 
m & Ai and 0, otherwise. 

We use lattice translation symmetry to block diago- 
nalize the Hamiltonian matrix The eigenstates 

that we use are the Bloch states^ 



1 



(2.2) 



1=0 



In this expression k is a reciprocal lattice vector (one 
of TV, where N is the number of sites), R; is a lattice 
vector (the order of I in the sum is not important), Ti 
is translation by R;, and A'nk is a normalization factor. 
The original basis states |n) are divided by translation 
into classes and any two states in the same class give 
the same Bloch state with an overall phase factor. A 
representative is chosen from each class and is used con- 
sistently to build the Bloch states. For a state \n) we 
denote its representative and the Hamiltonian ma- 
trix is block diagonal in the sense that (m, k'|_ff jn, k) = 
when k' 7^ k. The matrix (?fi, k|iJ|n, k) is diagonalized 
for each k vector using the well-known Lanczos method. 



B. State graph 

It is helpful conceptually to represent specific cases 
of our Hamiltonian by the state graph, in an abstract 
(or high-dimensional) space, such that each node corre- 
sponds to a basis state, and two nodes are joined by an 
edge of the graph if and only if t he two states differ by 
one particle hop. (As used in Sec. [HA, our many-body 
Hamiltonian is almost equivalent to a single particle hop- 
ping on the state graph.) The best-known precedent of 
this approach s Trugman's ptudy of one- hole and two-hole 
hopping in the t~Jz modelEH. 

Two states are connected if one state can be changed 
to another by a succession of hops. For states near half 
filling, our V = 00 constraint forbids many hop moves, so 
the graph is sparse and can possess interesting topological 
properties relevant to our spectrum. 



1. Droplets and accessibility of states 

We remarked (Fig. |l|) that an isolated hole is unable 
to hop; larger finite "droplets" of holes are still immo- 
bile, though they may gaiij kinetic energy from internal 
fiuctuations. Specifically,E3 if one draws a rectangle with 
edges at 45° to the lattice around the droplet, so as to 
enclose every site which deviates from the checkerboard 
order, then (with the V = 00 constraint) sites outside 
this rectangle can never be affected by fluctuations of 
the droplet. 

Due to the droplets, the state graph (representing all 
states with M particles) may be broken up into many dis- 
connected components, states which cannot access each 
other by allowed hops. This happens if both boundary 
vectors are even vectors, so that the system cell would 
support half filling with a perfect checkerboard order. 
For slightly less than half filling, the typical basis state 
consists of scattered immobile droplets of the kind just 
described. Each component corresponds to a different 
way to assign the holes to droplets, or merely a shifted 
way of placing the same droplets. j-Each sort of droplet 
has its characteristic energy levela23 and the system en- 
ergies are the sum, just like a system of noninteracting 
atoms, each having its independent excited levels. 

When the state graph is disconnected, the Hilbert 
space is correspondingly blocked into components which 
are not connected by matrix elements. This is beyond 
the blocking according to translational symmetry, which 
remains true. 



2. Building basis set from a starting configuration 

In systems with even-even boundary conditions, we are 
not mainly interested in the copfigurations with droplets, 
but in those with two stripes .c3 Since these are discon- 
nected from the rest of the Hilbert space, it speeds up the 
the diagonalization by a large factor if it is limited to the 
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stripe-configuration subblock of the Hamiltonian matrix. 
The only way to separate out this subblock is to gener- 
ate basis states iteratively from a starting configuration 
which has two stripes. 

What of the case of even-odd or (equivalently) odd- 
odd boundary conditions? These force one stripe (pos- 
sibly with additional holes, or possibly more stripes). 
In every such case (excluding L x L with odd L), the 
stripe can fluctuate and absorb any droplets, so we be- 
lieve the Hilbert space is fully connected, and the diag- 
onalization is not speeded up. Nevertheless, for every 
case near half-filling, we enumerated the basis from a 
starting configuration, as (we believe) this is the fastest 
way to do so. Systematically enumerating all possible 
fillings by brute force would take more time than the di- 
agonalization itself, since almost every brute-force trial 
configuration gets ruled out by the V = oo constraint. 

However, at low fillings (which we study in Ref. ||) 
the brute-force enumeration is used. In that case, if we 
generated from a starting configuration, there would be 
too many alternate routes to reach a given state. Hence 
much time would be wasted in searches to check whether 
each newly reached state was already on the basis state 
list. 

Below we outline the algorithm for building the basis 
states and the Hamiltonian matrix. 




"31 . 



FIG. 3: Building basis states from a given state n. Applying 
all allowed hoppings to n (the level state) gives the states 
ni, 712, M3, ..■ Applying hoppings to ni (a level 1 state) gives 
level 2 states nn, ni2, nia, ... The basis state list is then n, 
ni, n2, ns, nn, ni2, nia, n2i, 7122, nui, nii2, ... 
Note that when applying hoppings to a state, only add states 
not already on the list. We call this graph of states connected 
by hopping the state graph. 

To build basis states without using translation sym- 
metry from a starting state n, we follow the following 
steps. 

1. Apply hopping to n and obtain states ni, n2, n^,--- 
See Fig ^ Form the basis state list n, ni, n2, n^,... 
that are numbered 0, 1, 2, 3,... And record the 
Hamiltonian matrix elements Hq i, -ffo,2, ^o,3j--- 

2. Apply hoppings to the next state on the basis state 
list that has not been applied hopping to. Here 
from rii we get nu, ni2, ni^... Add these to the 
basis state list to form n, ni, 712, ris, nu, 7112, 



7ii3, ■•• (Note we should not add states that are al- 
ready on the list. For example, applying hoppings 
to ni will certainly give n back again. Do not in- 
clude this state in level 2 states. When computing 
Hermitian matrix elements we only need Hp q for 
p < q.) If 71 11 is the 7n-th element on the list, record 
the matrix elements Hi^m, Him+i, -ffi.m+2,--- 

3. Repeat the previous step until all the states on the 
list have been applied hopping to and no new states 
are created. 



"2, ns, 7111, 
And we have 



We finish with a basis state list n, ni 

ni2, ni3, 7121, "22, "111, ^1112, • 

stored the Hamiltonian matrix elements. 

When building basis states with translation symmetry, 
as in Eq. (2.2), the procedure is like the above, except 
we apply hopping to representative states of each class 
of translation related states and also store the represen- 
tatives of the resulting states. 



C. Boson and fermion statistics 

In this paper we study bo th th e boson and the fermion 
versions of our model Eq. (1.1), and one question that 
we often ask is when bosons and fermions have the same 
spectrum. In this section, we introduce a graph-based 
way to study the relationship between the boson and 
fermion spec tra. 

Using Eq. (2T), we know that if /3 is a state that can be 
obtained from the state a by one nearest-neighbor hop, 
then we have 



]H\a)f^sip.a){f3\H\a)b, 



(2.3) 



where the subscripts / and b denote the fermion and 
boson matrix elements respectively and we have used the 
fact that s{/3, a) = 1 for bosons. If we can write 



(2.4) 



i.e., s(/3, a) as a product of a function a that depends on 
one state only, then Eq. (2.3) can be written as a matrix 
equation. 



H/ — ^Ht)S, 



(2.5) 



where Sa,/3 — ^a,(30'{a). Because S — it is a simi- 

larity transformation, and then the eigenvalues of H/ are 
identical to that of Hft. Indeed, this is a much stronger 
condition than having identical eigenenergy spectra: the 
eigenstates are also identical (modulo a sign), so e.g. an 
operator which depends on the basis state has the same 
expectation in the boson and fermion cases. On the other 
hand, if the function a cannot be defined, the boson and 
fermion Hamiltonians are not equivalent. 

We can relate Eq. (p!4) to the loops which occur in 
virtually every state graph with the help of Fig. ^. Here 
s(/3, a) is defined on the arrow pointing from a to /?, 
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FIG. 4: Obtaining a{a), that is defined for each state, from 
a reference state and the step wise hopping sign function 
s{l3, a), that is defined on the arrow pointing from a to /3. De- 
fine cr(/3) — s{P,a)a{a). If the function (t(/3) is well-defined, 
i.e., all paths leading from to /3 give the same sign, then the 
boson and fermion spectra are identical. 



whereas the function <7{a) is defined on the node a of 
the graph. 

Fig. ^ suggests a natural way to construct a{a). First 
we choose a reference point. This should not matter and 
the one we choose is the starting state of the exact di- 
agonalization program, say 0. We set cr(0) — 1. Then if 
one hop takes a to /9, defin e (j{l3) = s(/?, a)a{a). This is 
certainly correct if Eq. (2.4) is true; but if the state graph 
has a loop, two different paths lead to the same state f3 
from 0, and cr(/3) can be well defined only if they produce 
the same sign. The answer is different for different state 
graphs (or disconnected subgraphs of the state graph), 
so it depends on the system dimensions and filling. 

We can check numerically whether the function (T{a) is 
well defined at the same time as we are building the basis 
set by applying hopping, i.e. constructing the state graph 
in the form of a tree (Fig. ||). We store a sign for each 
state, starting with a{n) = 1 for state n in Fig. ^. As we 
expand the tree, we compute the sign for the next level 
of states. Whenever we come to a state that is already 
on the list, we test whether the sign we produce following 
the current path equals that already stored for that state. 
Depending whether the sign always agrees or sometimes 
disagrees, we know the fermion and boson problems are 
or are not equivalent. This method is the basis for state- 
ment s late r in the paper on boson and fermion spectra 
(Sec. |III Al and Sec. |VI A| , and briefly in Sees. and 0); 
it is not necessary to compute all the eigenvalues. 

The boson-fermion equivalence condition can be re- 
stated in a way that is independent of the site order- 
ing convention at the start of Sec. Say each particle 
is originally numbered I, . . . , M and retains its number 
when it hops. Then (j{a, (3) is simply the sign of the 
permutation which changes the string of particles (when 
they are listed according to the site order) . The function 
a{a) is well-defined if and only if the product of cr(a,/3) 
around any closed loop of the state graph is unity; in 
other words, if and only if any sequence of hops that 
restores the original configuration induces an even per- 
mutation of the particles. When just two particles can 
be exchanged (an odd permutation), as is obviously pos- 
sible at low fillings, the fermion and boson spectra are 
obviously different. 



III. A SINGLE STRIPE 

It was first observed by Mila, in a model very similar to 
ours, that the creation of stripes is more favorable than 
of isolated, holes, as a way to dope the system below half 
filling.L3C3 We shall finally address more rigorously the 
phase stability of the stripe (Sec. VI] ), but in this section 



we will study the eigenstates of a single stripe in the 
system. In Fig. || we show a L^; = 4 and Ly — 7 system 
with 12 particles. A stripe of length four is formed along 
the X direction. The only four possible nearest-neighbor 
hops of the leftmost state are indicated with arrows and 
the resulting four states are shown on the right. Note 
that with a single stripe, the particles can only hop in 
the y direction. It is not possible for them to move in 
any fashion in the x direction. 
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FIG. 5; A single stripe in the 4x7 system. For the leftmost 
state, the four possible hops are labeled and the resulting 
configurations are shown on the right. The corre spond ing 
spin chain configurations are shown below (see Sec. IIIB). 



We also consider tilted boundaries (i^,5) x {0,Ly), 
which force a tilt of b steps to the stripe because the 
site (0,y) is identified with {Loc,y + b). In Fig. ^, the 
(5, 1) X (0, 7) and (6, 2) x (0, 7) systems are shown. Note 
because each stripe step is at 45 degree angle, an even 
(odd) horizontal length of the stripe must produce an 
even (odd) number of total vertical steps b. Thus L^ + b 
is always even. For a general {L^, b) x (0, Ly) system with 
Ly odd, the single-stripe states have M = Lx{Ly — l)/2 
particles. Note also that for the tilted stripes, vertical 
hopping is again the only allowed motion of the particles. 

Diagonal stripes are observed in La2Ni04j^ and they 
are an interesting topic in the t — J modeLca However, 
in our model with y = oo, no hopping is possible away 
from a 45° edge of any domain. Thus a diagonal stripe 
has no kinetic ener gy (un less there are additional holes, 
as discussed in Sec. |VD) and is disfavored. 



A. Boson and fermion statistics for one stripe 

We find from diagonalization that for one-stripe sys- 
tems with rectangular boundaries, the boson and fermion 
spectra are identical and are symmetric about zero. For 
a stripe in systems with tilted boundaries, there are com- 
plications. For the {2m -|- 1, 1) x (0, 21 + 1) system, with 
m and I integers, the spectrum is not symmetric about 
zero. In this case, with / odd, the fermion spectrum is 
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FIG. 6: A single stripe in the (5, 1) x (0, 7) and the (6, 2) x 
(0, 7) systems. The boundary conditions force a tilt of the 
stripe. The corresponding spin- chain s have total 5*2 equal to 
1/2 and 1 respectively (see Sec [II B). 



still identical to the boson spectrum, but for I even, the 
fermion energies are the boson ones vifith a minus sign 
(and because the spectrum is no longer symmetric about 
zero, the fermion energies are not the same as the boson 
energies). For the (2m, 2) x (0, 21 + 1) system, the boson 
and fermion spectra are identical and symmetric about 
zero. 

All these numerical findings can be e xplained in terms 
of the state graph, extending Sec. |ll C| . We must charac- 
terize the permutations induced by a sequence of particle 
hops which return to the same basis configuration. This 
is quite easy in the single stripe state, since (as already 
commented on Fig. ^ and Fig. 6|) , particles move in the y 
direction only. Our discussion includes all {Lx, b) x (0, Ly) 
systems (which are tilted boundaries if 6 7^ 0). In these 
cases, particles are confined to the same column, and the 
net permutation is a product of permutations in each 
column. 

Consider the total particle permutation induced when 
the system reLurns to the starting state by a loop on the 
state graph.cil If in this process the stripe as a whole 
makes no net movement across the y boundaries, then in 
each column the particles undo their hops and the net 
permutation is the identity. Thus, the only nontrivial 
loops are those in which the stripe crosses the y boundary. 
Indeed, since each hop moves the stripe locally by Ay = 
±2, the stripe must cross it must cross the boundary 
twice. The particles in each column have returned to 
the original positions with a cyclic permutation in each 
column, which has the sign 



(3.1) 



Therefore the sign of the system's permutation is crcoi^"" 
which (since L^ + bis even) reduces to (— l)''(^!'+i)/2: the 
fermion and boson spectra differ if and only if this phase 
is —1, i.e. only when b is odd and Ly = I (mod 4). 

The ±E symmetry of the single-stripe spectrum also 
falls out from visualizing the Hamiltonian as a single 
particle hopping on the state graph (with amplitude 
— t = — 1 for every graph edge), and recalling the spec- 
trum depends only on the attributes of closed loops. We 
can always change our convention for the phase factors 



of the basis states, which is equivalent to a "gauge trans- 
formation" on the nodes of the state graph. Whenever 
a state graph is bipartite, meaning every closed loop on 
it has an even number of edges, it is well-defined to di- 
vide the basis states into "even" and "odd" classes. If 
we change the phase factor by —1 for every "odd" basis 
state, every edge (every s(n, m) factor) picks up a factor 
— 1. By gauge equivalence, the new Hamiltonian matrix 
has the same spectrum. Yet on the other hand, it is 
manifestly the same as the old matrix with t —t, so it 
has the sign-reversed spectrum; this proves the ztE sym- 
metry. Now, the nontrivial loops in the state-graph are 
those that pass the stripe twice across the y boundary; 
this motion requires L^Ly = N hops. Thus, the state- 
graph is bipartite and the spectrum has ±E symmetry, 
if and only if N is even. 

If N is odd, the gauge-invariant effect of reversing 
t —t is to change the net product of s(n,rn) around 
every nontrivial loop by a factor —1. But we showed that 
switching Fermi and Bose statistics creates the same sign: 
hence the fermion spectrum is the inverse of the boson 
spectrum, as observed. 



B. Mapping to spin chain 

In Figs. ^ and ^, we have indicated a natural map 
from any configuration y{x) of the stripe to a state of 
a spin chain of len gth L^. (Recall that "configuration" 
was defined in Sec. [I A as a pattern of occupied sites; the 



map from the quantum states of the stripe is more subtle 
on account of phase factors, which are confronted later 
in this subsection.) Here it is obvious the spin length is 
one-half, for eaC|h-|Stripe step, y{x + 1) — y{x), can take 
only two values.eS The up step of the stripe is mapped 
to an up spin and the down step to a down spin. For 
example, the leftmost configuration in Fig. ^ maps to 
a spia state | itit)- This mapping was first noted by 
MilaEa, who used it (as we do here) to evaluate the exact 
ground state energy of a stripe. 

Note that when the periodic boundary conditions are 
rectangular, the stripe satisfies y{Lx) — y{0)- In the 
corresponding states of the spin chain, the total z spin 
component is zero. For tilted boundaries, the map to the 
spin-chain works exactly as before, the only difference 
being that the resulting spin configurations have — 
6/2. 

The stripe-to-spin-chain map is not one-to-one, for 
translation in the y direction of a single-stripe state gives 
the same spin-chain state (with a possible fermion sign). 
A rather analogous situation (but one dimension higher) 
appearS|_iij_the two-dimensional quantum dimer model 
(QDM).cZlo In the QDM, every dimer configuration cor- 
responded with a surface z(x,y), uniquely except for an 
arbitrary shift of z, just as our spin (or ID fermion) chain 
corresponds to a stripe y{x) here. There were two kinds 
of low energy excitations of the QDM: ripplon modes 
of the "surface" , and transverse motion (or tunneling) 
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of the "surface" through a periodic boundary condition. 
These correspond, respectively, to our ripplon modes la- 
beled by (see Sec. [II C 1 and also Sec. and to our 
y-translation modes labeled by ky (see Eq. ( 3.2S )). The 
analogy is imperfect in that, in the QDM, the dimer con- 
figuration is "real" and the surface is an abstract map- 
ping of it, whereas here the spin chain is the abstract 
mapping and the stripe is "real." 



mapping of the spin chain to a fermion chain (Sec. Ill C , 
below). In that representation, the phase factor can be 
associated with the hopping of a one-dimensional fermion 
across the periodic boundary, and it will be implemented 
by adding an artificial flux that pierces the ring for med b y 
the chain with periodic boundary conditions (Sec. [II D| ). 



Special case: zero gauge phase 



1. Spin Hamiltonian 

To determine the Hamiltonian in spin language, con- 
sider that each allowed particle hop takes the stripe either 
from up-step/down-step to down-step/up-step or from 
down-step/up-step to up-step/down-step. They corre- 
spond to the nearest-neighbor spin flips ti to and 
to ti respectively. 

Since the matrix element of H is — t between states 
differing by a hop, this must also be the matrix element 
of i?spin between states differing by a spin flip. Thus the 
equivalent spin Hamiltonian is 



spin 



i 

(-2i)5](^r5r+i+5r5,iVi)' (3-2) 



which is the so-called spin-1 
studied and exactly solvable. E 



'2 XX chain and is well- 



2. Mapping stripe states to spin states 

Lifted to act on Hilbert space, the map should take the 
mathematical form 



S\n) = 



(3.3) 



where \n) is a stripe state and jrig) is a spin state. 

Correspondingly, the equivalent spin Hamiltonian 
should satisfy 



SH\n) 



Hspin\Sn) 



(3.4) 



where H is t he sp inless fermion or hardcore boson Hamil- 
tonian, Eq. (1.1), and \n) is any single-stripe state. 

However, in the fermion case these equations are in- 
consistent with an arbitrary site-ordering scheme, for 
Eq. (3.3) did not take the sign factors into account. 
For many boundary conditions, there is no way to de- 
fine a unique function, on account of the large loops in 
the state graph illustrated in Fig. ^. That is, particles 
can be hopped successively (corresponding to a walk on 
the state graph) such that \n) returns to the state —\n). 
This gauge-like phase is familiar from the spinor coherent 
states (where each spin direction is mapped to a spin- 
1/2 ket), in which case a 27r rotation induces a —1 phase 
shift. This is handled most transparently after a further 



The many-to-one nature of the map (because j/ — > j/+2 
translation of the stripe produces the same spin state) is 
easier to handle. Let us choose a rectan gula r system, 
adopt the site-ordering conventions of Sec. II A, and pro- 
visionally assume there is no gauge-like phase ambiguity. 
We ask that if Tr translates a state in the y direction by 
Ry, then STu^ \n) = S\n). One way to recover the notion 
of a unique map is to form Fourier states like Eq. ( ^.2| ), 
but only in the y direction: 



\n,ky) = 



(3.5) 



Ly because 



The normalization factor here is always 1/ 
the Ly translations in the sum cannot produce two identi 
cal states. From group theory we know (m, k'y\H\n, ky) = 
ii ky ky, so states with different ky belong to discon- 
nected blocks. In each block with a given ky, there is a 
one-to-one correspondence \n,ky) ^ S\n). 

Using Eq. (2.1), because H commutes with translation 
Tr , it is easy to show that 



H\n,ky)~{—t) s{n,m)\m,ky 



For q e M, we then have 



(p, ky\H\n,ky) = -s{n,p)t. 



(3.6) 



(3.7) 



But if we translate p by Ry so we have Tn,^\p) = cr\p'), 
where a is the fermion sign from translating the state p 



by Ry, then the translation state \p' , ky] 
Then we have. 



{p',ky\H\n,ky) = -s(n,p)(Te-''=«^n 



\P, ky) 



(3.8) 



Eqs. ( |3.8[ ) and ( p.7| ) tell us that choosing p' rather 
than p, the matrix element {p' ,ky\H\n,ky) has an ex- 
tra sign a and an extra phase factor e^'*^"^". However, 
in the spin-chain language, we have, {Sp\Hspin\Sn) — 
{Sp' \Hspin\Sn) , because S\p) — S\p'). This means that 
we should have ae~^''^^ 
(mod 4), then, as we have shown in Sec. 
fermion hop in the one-stripe case produces a -1-1 sign,Eil 
thus we have s(n,p) = 1, and at the same time, every y 
translation gives cr = 1. In addition, if we choose ky ^ 0, 
then we have, 



1 for all Ry. Choose Ly 



HI A, evee 



{p,0\H\n,0) ^ -t. 



(3.9) 
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for any \p) and |n) in their respective translation classes. 

In summary, for ky = and Ly = 3 (mod 4), we 
map the Fourier transformed single-stripe state \n, 0) to 
the spin-chain state \Sn), then we have an identity of 
matr ix e lements {p,0\H\n,0) = {Sp\Hspin\Sn) , and thus 
Eq. (3^) is satisfied. 

For general ky 7^ or Lj, ^ 3 (mod 4), this map in- 
volves additional phase factors which will be discussed 
in Sec. HID using a further map to one-dimensional 



fermions that wc will describe in Sec. [II C 



We have checked that the eigenvalues and eigenvec- 
tors of the spin chain system match those of the spinless 
fermion and hardcore boson problems, using a spinjiiag- 
onalization program developed for another projecto 



C. Fermion representation of a stripe 

It is possible, in turn, to map each configuration of 
"spins" to that of a one-dimensional (ID) lattice gas of 
particles: every | becomes an occupied site, every | be- 
comes a vacant one. The total number of "up" steps, 
i.e. fermions, is = [L^ + 6)/2 if the boundary vector 
along the stripe is (L^,, 5). The +y hop of a real particle 
on the square lattice (implying a —y fluctuation of the 
stripe), translates to a spin exchange ti^it in the XX 
chain, and finally tQ_a +x hop of a ID particle. 

As is well-knownjlj in one dimension hardcore parti- 
cles may always be treated as fermions: if there is no 
path for them to exchange, then the statistics has no 
physical meaning. One can describe the system as one- 
dimensional spinless fermions which are noninteracting, 
as the Fermi statistics is already sufficient to keep two 
particles from occupying the same site. The explicit re- 
lation of the stripe path y[i) to the ID spin and fermion 
representations is 

yii + l)-yii) = 2S,, = 2n,-l. (3.10) 

where hi = c[ci. The effective one-dimensional Hamilto- 
nian is 

'n=-tY,{4+lC^^+4c,+^), (3.11) 
i 

thus the dispersion is 

-2tcos9, (3.12) 

where q is the one-dimensional wavevector. We can con- 
struct exactly the ground state by filling the lowest- 
energy plane-wave states, up to the one-dimensional 
Fermi vector vt: {— 7r/2 for an untilted stripe); here 
V — Nj^/Lx is the density of one-dimensional fermions, 
corresponding to a stripe slope 2v —1. All excited states 
correspond exactly to other ways of occupying the ID 
fermion states. 



1. Excitations of a stripe 

If a stripe is coarse-grained as a quantum-fluctuating 
string, then it can be approximated using a harmonic 
Hamiltonian, 



n 



— dx 



E 
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dy 



ax 



p2 

g 



KqWa 



(3.13) 



Here p* is the effective mass density per unit length of 
the stripe, and the stiffness K is analogous to a string 
tension. The long-wave excitations of such a "string" are 
the quantized capillary waves known as ripplons, which 
have dispersion 



c\q\ 



with 



(3.14) 



and q is the one-dimensional wavevector. 

On the other hand, in the microscopic representation 
by ID-fermions, the fundamental excitation is evidently 
a particle or hole, which corresponds to a mobile kink of 
the stripe. The ripplon is thus a composite excitation, a 
kink-antikink bound state. In the equivalent language of 
the spin-1/2 XX model, the ripplon maps to a magnon 
(hydrodynamic mode), while the kink maps to a spinon. 
The fractionalization of the spin-1 magnon into two spin- 
1/2 spinous is a familiar fiact, as this system is a special 
case of a Luttinger liquid. c£l 

We can use the F ermi sea representation to extract the 
parameters in Eq. ( 3.13 ). As noted above, the ID Fermi 
wavevector is (1 -|- dy/dx)TT/2 and by Eq. ( 3.12| ) the total 
energy is ^^^{—2t)cosq, where this notation means 
the sum over occupied ID fermion states. Equating the 
energy density due to small tilts to ^K{dy/dx)'^, we ob- 
tain K ~ TTt/2. Next, the ripplon velocity is the same 
velocity as the Fermi velocity, c ^ vp — 2t/h. (This 
is a stand ard fact of onc-dimensional Fermi seas.) With 
Eq. (|3.14D, that implies 



2. Stripe roughness 



(3.15) 



The ID Fermi sea representation also allows the exact 
computation of stripe fluctuations as a function of x. It is 
straightforward to show (after rewriting in terms of wave 
operators) that at half filling, 

1, ifr = 0; 

{{2h, - l){2h,+r - 1)) = <( 0, for even r ^ 0; 



for odd 



We also know from Eq. ( 3.10 ) that 

R-l 



(3.16) 



y{i + R)-y{i) = J2i^h,-l). (3.17) 



3=i- 
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This and Eq. (3.16) give that 



R-l 



TT^ ^ — ' \ r 

r odd— 1 



1 R 



const H — ^ In i? 



(3.18) 



Not surprisingly, the same result can be derived via the 
zero- point intensities of the harmonic ripplon modes, of 
Eqs. (ETl and (|]l|). 

Eq. (3.18) mea ns that, due to the anticorrelations evi- 
dent in Eq. ( |3.16 ), the transverse deviations of the stripe 
grow slowly with length. Thus, if it moves transversely 
over larger distances, the stripe can often be approxi- 
mated as a single rigid object. This approximation will 
be invoked in several sections to explain how energy split- 
tings depend on the system size Ly in the transverse di- 
rection. 



D. Phase factors and stripe effective mass 



The mappings described in Sec. Ill C| are less trivial 
than they appeared, due to two related facts. The first 
fact is that the mapping of stripe to spins (or to ID 
particles) is many-to-one, as already noted in Sec. Ill E . 

The second fact is that the statistics of ID hard- 
core particles is not quite as irrelevant as suggested in 
Sec. Ill C: they can be permuted by moving them in the 
X direction through the periodic boundary conditions. 
As is well known, this induces additional phase factors 
in a finite system. Consider, for example, a sequence of 
—y hops along a stripe, such that exactly one particle 
hops in each column. The net effect on the real particle 
configuration is to translate the stripe by (0,2); in the 
fermion case, a sign factor ct^qj is also picked up, where h 
is the number of columns in which t he r eal particle hop 
crossed the cell boundary, and Eq. (3^) is the fermion 
sign picked up from the resulting rearrangement of cre- 
ation operators among sites in a column. (Note that, 
here and for the rest of this discussion, we limit consider- 
ation to lattices in which the second boundary vector is 
(0, Ly).) Meanwhile, the same sequence of real particle 
hops maps to a cyclic permutation of the ID particles in 
the -\-x direction - exactly one -\-x hop occurs on every 
bond along the chain. 

If the wavevector of our eigenstate is taken as (fc^^, ky)^ 
the stripe configuration shifted by (0, 2) ought to have an 
amplitude in the wavefunction e*^'"'" times the amplitude 
of the original configuration. In a ID fermion system, 
however, the actual phase factor is (—1)^+"^, i.e. the 
number of fermions that are cyclically permuted. We 
can account for all the phase factors by modifying the 
one-dimensional chain so that a particle hopping around 
its periodic boundary condition picks up a phase e*"^ in 



other words by inserting a flux into the ringjlj where 
(p = 2ky + (iV+ — l)7r is sufficient to account for the two 
phases mentioned above, in the boson case; in the fermion 



case, the an additional term jLx{Ly + 1)it is needed.i^l 
Finally it is equivalent, modulo 2-k, to replace by b in 
the last term: 



1, 



= 2ky + (iV+ - l)7r + [-b{Ly + l)7r] 



fermions •, 



(3.19) 



where the last term is included only in the fermion case, 
and even then it is zero in a rectangular system. [Notice 
that Ly is always odd when the boundary vector is (0, Ly) 
transverse to the stripe, and the other boundary vector 
{Lx,b) must be even.] It follows that the allowed ID 
wavevectors are (27rm -I- (f)/ L^, where m is any integer. 

We get = when the number of ID fermions 
is odd and <j) = -k when is even; either way, the ID 
fermions occupied in the ground state are placed symmet- 
rically about q = 0. In any eigenstate, the real wavevec- 
tor component is = occ?' ^^vs k^ = Q for the stripe 
ground state in any rectangular system. It is not hard 
to compute occ(^^^) cosg to obtain the (exact) total 
ground state energy 



£^stripc(0,0) = -2Vsin(Vi,), 



(3.20) 



which implies the stripe energy is —2t/pi per unit length 
in the thermodynamic limit All ou r numerical results 
for stripe states agree with E q. (|3.2(]| ) (see Table |) . 

When, as mentioned in Sec. Ill C 1| , we model the stripe 
direction as a free particle moving in one dimension, it 
is useful to know its effective mass m* . The minimum- 
energy state of a single stripe with ky ^ Q consists of 
rigid motion in that direction, so for small ky we expect 



n^k 2 

-£'stripc(0, A;„) — £'stripc(0, 0) = — — . 

2m* 



(3.21) 



In the ID fermion representation, this state is con- 
structed by shifting the ground state Fermi sea by a 
wavevector 2ky/Lx^ a s foll ows from the first term in 
the phase shift Eq. ( 3.19 ). From this one obtains 
L^stripc(0, fc „) = £'striDc (0,0) cos{2ky/Lx). Combining that 



with Eqs. ( 3.20 ) and ( 3.21 ), one obtains 



2t2 



m 



h'L\ 



■ sin(7r/ij;). 



(3.22) 



For a long stripe, Eq. ( |3.22D implies an effective ma ss de n- 
sity p* = limL^ocTT^* /Lx in agreement with Eq. (3.15). 

To obtain the effective mass m* from exact diagonal- 
ization data, we use the smallest ky = 2t: /Ly with a 



large Ly = 101. Eq. ( 3.21 ) gives the following expression 



to extract m* numerically. 



{2n/Ly)' 



2{Estripe{0, 2'K I Ly) — EstripeiO, 0)) ' 



(3.23) 



The resul ts ar e shown in Table |[ they agree perfectly 
with Eq. ( |3.22D . 

It should be noted that the above calculation of the 
phas es , is entirely equivalent to that of Green and Cha- 
monOtJ (where it is applied to a stripe on a triangular 
lattice), except that we include the case that the particles 
are fermions. 
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TABLE I: Stripe effective mass m* [Eq. (3.23)] using single- 
stripe excited state energies with k — (0, ky), with Ly — 101. 





Estripe{0,0) Estripe{0,2Tv/Wl) 


m* 


4 


-2.8284271 


-2.8270590 


1.4143276 


6 


-4.0000000 


-3.9991400 


2.2500806 


8 


-5.2262519 


-5.2256198 


3.0615292 


10 


-6.4721360 


-6.4716350 


3.8627623 


12 


-7.7274066 


-7.7269913 


4.6587845 


14 


-8.9879184 


-8.9875635 


5.4517988 



hole-with-a-stripe spectra are identical. The numerical 
method described in Sec. II C is used to check whether 



the function a{a) is well defined using the state graph 
Fig. For a Lx X Ly system with even and Ly odd, 
the one-hole-with-a-stripe state has M = {Lx — l)Ly/2—l 
particles. For the one-hole systems we checked, for ex- 
ample, 4x9 with M = 15, 4 X 11 with M = 19, and 6x7 
with M = 17, a{a) is always well defined£d 



A. Energy dependence on Ly 



IV. ONE HOLE ON A STRIPE 

For a system with a single stripe, we have seen that 
the particle motion is strictly limited to sliding in the 
direction perpendicular to the stripe, and this enables us 
to map our two-dimensional system to a one-dimensional 
spin chain and solve it exactly. With holes added to the 
stripe state, a number of new motions are allowed. In 
Fig. ^ we show a 6 x 7 system with a stripe and a hole. 
We see that the hole can now move along the stripe and 
also, the stripe can fluctuate and leave the hole stranded 
(i.e., immobile). In this section we study the onc-hole- 
with-a-stripe problem. 
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FIG. 7: A 6 x 7 system with a stripe and a hole. The hole is 
indicated in the box. It is initially lying on the stripe (leftmost 
state). The arrows indicate the particles hops that result in 
the hole moving along the stripe to the right (middle state). 
It is also possible for the stripe to fluctuate and leave the hole 
immobile (rightmost state). 

Motion of holes has been analyzed on domain walls 
in the anisotropic t-J model (the t-Jz model with sorae 
added terms), using essentially analytical techniques. E2l 
The most interesting observation in Ref. ^ is that the 
charge-carrying "holon" is associated with a kink of the 
wall, and thus carries a transverse "flavor" (reminiscent 
of our model's behavior in Sec. IV D). A hole on our 



stripe differs in an importsnt way from the holes on a 
stripe in the t-J^ model.ESE] First, in our case the refer- 
ence stripe alre ady has 1/2 hole per unit length, whereas 
in Refs. ^ and |30| it is a plain domain wall without holes. 
Second, Jz/t was not too large in those models, so that 
hopping of the hole several steps away from the stripe 
has a noticeable contribution in the wavefunction, while 
it is unimportant in our case. 

^From diagonalization, we observe that for systems 
with rectangular boundaries, the boson and fermion one- 



When a hole is added to a stripe, more hops are allowed 
and the state gains kinetic energy, so the energy is lower 
than that of the single stripe of the same length. We 
define the energy difference. 



^{Lix, Ly) = Ehole{Lx, Ly) — Estripe{Lx), 



(4.1) 



where EhoieiLx, Ly) is the ground state energy of one 
hole with a stripe on a L^; x Ly lattice and Estripe{Lx) 
the ground state energy of a single stripe with length Lx- 
Here E^oieiLx, Ly) is the same for bosons and fermions, 
as we discussed abov e, and Estripe{Lx) does not depend 



on Ly (see Eq. (^^). 

A plot of A{Lx,Ly) vs Ly shows a fast decay in Ly so 
we try the following exponential fitting function. 



A{Lx,Ly) = K{Lx)-A{Lx 



(4.2) 



where h{Lx) = A{Lx,oo), A{Lx), and l{Lx) are fitting 
parameters that depend on the length of the stripe Lx- 
(We wi n inve stigate the dependence on Lx of A{Lx) later 
m Sec. |IV C| .) We choose a minus sign in front of A{Lx) 
because A{Lx,Ly) < A{Lx)] in Eq. (4^), A{Lx) is posi- 
tive. 



This fitting form Eq. (O), suggests the following linear 
regression check. 



ln(A(L,, Ly + 2) - A{Lx, Ly)) = C - 



liLxY 



(4.3) 



where C is a constant that depends on A^ I, and Lx, but 
not on Ly. In Fig. ^ we plot ln(A(La;, Ly+2)~ A{Lx, Ly)) 
vs Ly for Lx = 4, 6, 8, 10, 12, 14 and Ly = 5, 7, 9, .... The 
linear fit is excellent for all data sets. 



B. Stripe potential well and effective mass 

Ho w c an we account for the exponential fitting form, 
Eq. (^^)? Let us consider the hole fixed at some posi- 
tion and the stripe meandering in y direction. In Fig. |^, 
we have shown that the hole can be in contact with the 
stripe or the stripe can fluctuate away, leaving the hole 
behind and immobile. When the stripe is in contact with 
the hole, the energy is lower than the energy of a single 
stripe Estripe, which is also the energy when the stripe 
is separated from the hole. Because we have periodic 
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''stri ve , 



The slope for each Lx curve is —l/l{Lx) in Eq. (1.3) 



boundary conditions in the y direction, we can use a pe- 
riodic array of potential wells to model the y motion of 
the stripe; the well is at the position of the hole. 



symmetric in one well). 

As far as our hole-with-a-stripe problem is concerned, 
the well separation d is Ly, the barrier height Vb is 
E stripes the ground state energy of one isolated well i?o is 
Ehoie{Lx, Ly = oo), and the ground state E2 of the sys - 
tem with tunnehng is Ehoie{Lx, Ly). Therefore, Eq. (4.4) 
for the wells with tunneling translates into the following 
equation for our hole-with-a-stripe problem, 

Ehoie{Lx,Ly) - EhoieiL^,^) = -A(L,)e-^«/'(^==). 

__ (4.5) 

Using th e de finition for A{Lx,Ly), Eq. ([4.l|) , we see 
that Eq . (4.5) is exactly the fitting form we used before, 
Eq. ( |4.2|) . In addition, the inter- well tunneling amplitude 
Eq. (Wm gives us a way to calculate the effective mass 
m*{Lx) of the stripe of length L^, 



j—r = \J 2m* {Lx){Estripe{Lx) - Ehoie{Lx, 00)), (4.6) 



Eq. (4.1). where we have set S. = 1. We get 



m*{Lx 



l/l\Lx) 



2{Estripe{Lx) — Ehole{Lx,Oo)) 



(4.7) 



Using the linear fitting slopes in Fig. |8| (tha t are 
— l/l{Lx)), we can compute m*{Lx) using Eq. (|4.7|), and 
our results are shown in Table The effective mass 
results are consistent with that obtained from the single- 
stripe energy dispersion relation in Table ||. 




FIG. 9: Potential- well and barrier potential used to un- 
derstand the exponential decay of one-hole energy in Ly in 



Eq. (4.2). y is the stripe position in the direction perpendic- 
ular to the stripe, Vb is the depth of the well, Eq the ground 
state energy of one isolated well, E2 < Eo the ground state en- 
ergy including the tunneling between wells (symmetric state), 
and d the separation betw een two wells. E2 — Eq decays ex- 
ponentially in d as in Eq. (4.4). 



In Fig. H we show potential wells separated by a bar- 
rier of thickness d. Eq is the ground state energy of 
an isolated well, and E2 the ground state energy of the 
multiple-well system, which corresponds to a symmetric 
state and is lower than Eq- From standard quantum me- 
chanics textbooks (see e.g., Ref. ^2|), we know the ampli- 
tude to tunnel between adjacent wells difference decays 
exponentially with their separation. 



E2-Eo = -Ae-'^l\ where, 



V2m(Uo - E^) 



(4.4) 

where m is the mass of the particle in the well, and the 
constant A is positive (because the wavefunction of Eq is 



TABLE IL Stripe effective mass m* calculated from Eq. (4.7) 
using the potential-well model for the one-hole-with-a-stripe 
problem. liLx) is the decay length in Eq. (4.2) obtained from 
linear fitting in Fig. ^ Ehoie{Lx, Ly) is used to approxmate 
EhoieiLx, 00) using the large Ly listed in the table (as the ex - 
ponential decay of Ehoie{Lx, Ly) in Ly is fast, see Sec. IV A[ ); 
and m*{Lx) is the effective mass of the length-Lo, stripe cal- 
culated using Eq. (1/7). Estripe{Lx) has appeared in Table ^ 



Lx Ly Ehoie{Lx, Ly) l{Lx) m*{Lx) 



4 25 
6 25 
8 25 
10 21 
12 19 
14 17 



-3.00000009 
-4.29850460 
-5.60276427 
-6.89914985 
-8.18958713 
-9.47602246 



1.4659 
0.8724 
0.6648 
0.5553 
0.4868 
0.4388 



1.3561 
2.2009 
3.0050 
3.7979 
4.5652 
5.3197 



The main physical significance of the result in this sec- 
tion is its implications-|£oi' the anistropic transport of 
holes in a stripe arrayflCj The conductivity transverse 
to the stripes depends completely on mechanisms which 
transfer holes from one stripe to the next. The high en- 
ergy of the "stranded" state, in which the hole is immo- 
bile and off the stripe, means that holes are mainly trans- 
ferred when the stripes collide, either directly at contact 
or else by a delayed transfer, such that the hole spends 
a short time in a stranded state until the second stripe 
fluctuates to absorb it. 
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C. Hole dispersion and Lx dependence 



In this section we fix Ly and study the dependence of 
A{Lx, Ly) on the length of the stripe L^. In Fig. |l^ we 
plot the energy difference A{Lx,Ly) for {Lx,b) x (0,7) 
systems with 6 = 0, 1, 2.eJ 
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FIG. 10: A{Lj:, Ly) = Ehoie — Estripe as a function of l/L^ 
for three classes of lattices {Lx,b) x (0, 7) with fe = 0, 1, 2. 

We see fro m Fig . |l0| the following fitting function (jus- 
tified in Sec. IVD) works well for both rectangular and 
tilted lattices, 



A(L,) 



C{b) 



(4.8) 



This fitting form enables us to extrapolate the energy 
gap formed by adding one hole to an infinitely long (Lx) 
stripe. The intercepts of the three curves for the three dif- 
ferent classes of lattices all approach A = —0.66. Later, 
in Sec. VII B , we will study the stability of an array of 
stripes and we will use the value for A because we want to 
know whether holes added to a stripe stick to the stripe 
to form a wide stripe or a new stripe. A is energy lowered 
by adding a hole to a stripe and will be relevant there. 
It is a binding energy in the sense that one hole off the 
stripe (immobile) has zero energy. 

A hole is mobile along a stripe; one expects its disper- 
sion relation to be 



2ml 



(4.9) 



on a long stripe oriented in the x direction, where q is the 
hole's wavevector. We can estimate the effective mass in 



Eq. (4.9) from the numerics, if the lowest-energy state 
with system wavevector k = (q, 0) is produced by boost- 
ing the hole to wavevector q. At small g, this has a much 



smaller excitation energy than the ripplons (stripe ex- 
citations) which have linear dispersion [Eq. ( 3.14 )]. In- 
deed, the numerical spectra for systems with L^ — 10 
and Ly — h show these contrasting dispersions for a 
stripe with and without a hole. Fitting the difference 
E{q) - E{0) to h^/ml{l ~ cosg), as a plausible guess, 
we estimate fi} I2m\ « 0.3i; the effective mass m\ of a 
hole bound to a stripe is thus about six times as big as 
that of a free particle in an empty background (for which 
r^llm* = 2t). 



D. Hole and stripe steps 

How does the hole interact with the stripe fluctuations? 
In some models, hole hopping is suppressed if the stripe 
is tilted away from its favored direction. In that case, the 
"garden hose effect" is realizediES stripe fluctuations are 
suppressed when the hole is present, so as-^o optimize 
the kinetic energy of the hole on the stripeE3 

However, in our model, a stripe tilt actually enhances 
hole hopping. This is clear in the extreme case of a stripe 
at slope +1, as occurs at maximum filling in the L x 
L square system with L odd: the bare stripe has zero 
energy. When the system is doped by placing one hole 
along an edge of the stripe (see Fig. |ll), the accessible 
configurations are equivalent to a one-dimensional chain 
of the same model, with length 2L and L — 1 particles 
on it. The one-dimensional chain has two domain walls, 
one of which appears as a down-step of the stripe, the 
other of which appears as an up-step bound to a hole on 
the stripe. In the limit of large L, each domain wall has 
kinetic energy — 2t. The total energy is — 4t, or about 
seven times larger than the energy difference —A for a 
stripe with no tilt. 

This picture is supported by the trend of the stripe en- 
ergy with one hole, in the presence of a boundary vector 
{Lx,b) that forces a tilt. For L^, = 5 to 8, the minimum 
energy, as a function of b, occurs around 2.5. [The mea- 
sured minimum is at 5 = 2 in the case that L^ and b are 
even, at 6 = 3 in the odd case, while the second-lowest 
energies are at 6 = 4 and 6=1, respectively] This is 
consistent with the_idea that the hole on a stripe binds 
with some up stepsEll which total boi where b^ 2.5. 

This would imply that holes tend to come in "up" and 
"down" flavors; this is not an actual quantum number, 
since the "steps" conside red here (in contrast to "kinks" 
mentioned in Sec. Ill C 1 ) are not quantized objects. Our 
hole-step binding has a similar origin as that in Ref. |2^: 
the local configuration is the same as a fragment of stripe 
rotated 90° from the main stripe. 

A corollary of the hole-step binding is that a hole 
should force an extra tilt \b — bQ\/ Lx on the rem aining 
portion of the stripe. Thus the tilt energy in Eq. ( 3.15 ) 
contributes a size dependence C(b)/Lx, where 



1 



C{b) = -K[{b - bo)' - b'] = Kboibo/2 - 6), (4.10) 
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FIG. 11: A 7 X 7 system is shown, at one particle less than 
maximum filling. All motions in this system occur within 
the strip bounded by dashed lines, which is equivalent to a 
14 X 1 chain of spinless V = oo fermions. The stripe edge is 
indicated by thick solid lines. 



that is 4.9, 1.0, and -2.9 for 6 = 0, 1, and 2 respectively. 
This is consistent with the clear l/L^ size dependence 
we found numerically (see Fig. |l^); the fitted coefficients 
for 5=1 and 2 are in re asonable agreement with C(6) 
as predicted by Eq. ( 4.10 ). (Presumably C(0) is reduced 
from our prediction because the hole is not always, or not 
usually, bound to an up-step in the case b — 0.) 



V. TWO HOLES ON A STRIPE 

When two holes meet on a stripe, they create enough 
room that particles can exchange with each other, so the 
boson and fermion energy spectra are no longer the same. 
As in the one-hole case, we study the energy of the two- 
holes-with-a-stripe problem as a function of the two di- 
rections of the lattice. 



A. Energy dependence on Ly 



As in Eq. ( fl.lj) for the one-hole case, we define for the 
Lx X Ly lattice the energy difference between the case of 
two holes with a stripe and that of a single stripe, 



A2{Lx,Ly) = E, 



2hole 



where the subscript 2 denotes the two-hole case. 
(Strictly speaking, we should write /\^2^ {Lx,Ly) and 

^■z'hoiesi^x, Ly) because these quantities are not the same 
for boson and fermion cases. Here, without the super- 
scripts, they stand for both cases.) Plots of A2{Lx, Ly) 
vs Ly for bosons and fermions show fast decay similar 



AL, 



(5.1) 



to that in the one- hole problem in Sec. IV, except here 
the boson and fermion curves approach different values 
at large Ly. 



As in Eq. (4^) for the one- hole problem, we write 



A2(ix, Ly) = A2(L^) - A2iLx)e' 



(5.2) 



\n{A2{Lx, Ly 4- 2) — A2{Lx, Ly)) vs Ly for bosons and 
fermions have been plotted with L^ = 4,6,8,10. The 
exponential dependence is checked nicely as in Fig. ^ for 
the one-hole case, and the slopes are obtained and will be 
used for calculating the effective mass in the next section. 



B. Stripe potential well and effective mass 



The resemblance of our treatment of the two-holes- 



with-a-stripe problem in Sec. VA with that of the one- 



hole problem in Sec. IV A prompts us to ask whether 
the two-hole problem can be studied using a one- 
dimensional e ffecti ve potential like the potential wells 
used in Sec. 



IV B 



Here with two holes, we have a 
more complicated problem because the relative positions 
of the two holes and the stripe can have three cases: 
two holes on the stripe (with energy E2hoies{Lx, oo)), 
one hole on the stripe with one isolated hole (with en- 
ergy Eihoie{Lx, oo)), or one stripe with two isolated holes 
(with energy E^tripeiLx))- 

As an approximation, we consider the deepest well 
{E2hoies{Lx, oo)) Only because that is the most proba- 
ble position to find the particle. Then the effective mass 
equation Eq. ([4.7[) is modified to become. 



m*2{Lx) = 



l/lliLx) 



'^{^^ stripe {Lx) — E2holes{Lx, oo)) 



(5.3) 



where hiLx) comes from the linear fitting slopes. The re- 
sults for from this two-hole-with-a-stripe calculation 
are in Table pi^^ They are comparable to the one-hole 
results in Table ||. 



TABLE III: Stripe effective mass 



calculated from 
Eq. (5. J) using the double- well potential model for the two- 
hole s-wit h-a-stripe problem. hiLx) is the decay length in 
Eq. (5.2) obtained from linear fitting. E2hoies{Lx, Ly) is used 
to approxmate E2hoies{Lx, oo) using the large Ly listed in the 
table; and mKix) is t he e ffective mass of the length-L^, stripe 
calculated using Eq. (5.3). The superscripts b and / denote 
bosons and fermions respectively. Estripe{Lx) has appeared 
in Table ^ 



2holes 



E. 
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1X12 



7712 



4 19 -3.896952 -3.818556 0.5919 0.6135 1.3358 1.3417 

6 19 -5.238613 -5.204899 0.4317 0.4386 2.1664 2.1574 

8 19 -6.524636 -6.510024 0.3621 0.3643 2.9375 2.9349 

10 15 -7.797992 -7.790359 0.3167 0.3274 3.7605 3.5377 



Energy dependence on Lx 
binding 



and stripe-step 



To compare with Fig. |10| for one hole with a stripe, in 
Fig. |l^ we plot A2{Lx, Ly), for three classes of lattices 
(Lx^b) X {0,Ly = 7) with b = 0,1,2, vs l/Lx- Here the 
1/Lx 



fit, as in Eq. (1 



for the one-hole problem, is no 
longer good. But we can still extrapolate the energy gap 
of the two-hole state and the stripe state for an infinitely 
long stripe. For the {Lx,0) x (0,7) lattices, the gap is 
extrapolated to be —1.4 approximately, for bosons and 
fermions. Comparing to the single hole gap of —0.65, 
we see that binding between holes are not strong. This 
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information will be useful when we consider stripe-array 
formation. 
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FIG. 12: A2(L. 



Ly) — E2holes 



-E stripe as a function oil/ 



for three classes of lattices {Lx, b) x (0, 7) with 6 = 0, 1, 2. 

The size dependence of two holes on one stripe can 
be interpreted by the notion (Sec. iVD) that holes are 
bound to up- or down-steps. When the boundary vec- 
tor is {Lx,b) with 6 = or 1, the stiffness cost will be 
minimized if the holes adopt canceling flavors "up" and 
"down" (so long as 5 < 2bo « 5). This may explain why 
the 6 = curves in Fig. [ij seem to have a weak coeffi- 
cient of l/Lx. The fact that A2 « 2Ai indicates weak 
interaction between the holes, and that can only happen 
if the holes (of opposite flavors) tend to repel. 

However, on reflection it will be noticed that, in a suf- 
ficiently large system, added holes can be "condensed" 
to create segments of stripe oriented in the ±y direction. 
These "90°" segments will surely appear in the ground 
state of one stripe in a large system, as the energy per 
hole to increase the net stripe length is far lower than the 
energy to add a free hole hopping along the stripe. Thus, 
holes do attract, on a sufRciently long stripe. That does 
not contradict our numerical observation of repulsion: as 
just explained, the constraint of zero net kink "flavor" on 
a stripe forces two holes to adopt opposite flavors, which 
repel, whereas the holes that form a single 90° segment 
all have the same flavor. 

Two effects compete with the formation of 90° seg- 
ments in the case of short stripes or small numbers of 
added holes: (i) Say the stripe contains one 90° segment 
with h holes; that extends ±2h {not ±boh) in the y direc- 
tion, which forces the remainder of the stripe to have a 
slope ^2h/Lx and the associated tilt energy Eq. ( 3.13 ). 
(ii) Say the stripe contains two 90° segments, with the 
opposite directions: no net slope is forced, but twice 
as many 90° corners are present, and presumably each 
90° bend costs energy (since it suppresses stripe fluctua- 
tions). 



The implication is that the 6 = and 6=1 curves 
in Fig. |l^ are heading towards a well-defined asymptote 
at — 2Ai, corresponding to a pair of repelling, opposite- 
flavor holes. However, at sufficiently long L^, they must 
cross over to a different curve with a lower asymptote, 
corresponding to a pair of attracting, same-flavor holes, 
bound into an incipient 90° segment. We conjecture that 
the entire 6 = 2 curve is in the latter regime, while the 
downturn of the 6=1 curve at the largest suggests it 
is beginning to make the crossover. 



VI. TWO AND MORE STRIPES 



Starting from this section, we study the interaction 
among stripes. We use the same diagonalization pro- 
gram for the problem of one stripe and one stripe with 
holes that was introduced in Sec. ||. Again we do not 
exhaustively enumerate all possible states with a given 
number of particles M. Instead, we build basis states 
from a starting state which has the stripes all merged 
together. The motivation of these studies is that under- 
standing of the stripe-stripe repulsion is a prerequisite to 
calculating the stripe-array energy as a function of filling 
near n = 1/2, which in turn is needed in order to ascer- 
tain the density at which sthe stripe array would coexist 
with a liquid phase. Unfortunately, we can only study 
comparatively short stripes and the asymptotic form for 
La; ^ 00 is quite unsure from our numerical results. 



A. Boson and fermion statistics 

We observe from exact diagonalization that for rect- 
angular lattices withJavo stripes, the boson and fermion 
spectra are identicaLEa In Fig. |l^, we show two configu- 
rations of two merged stripes. For a particle on the stripe 
we see that the horizontal hops are limited by nearest- 
neighbor replusion to adjacent columns only. This means 
that, with two stripes, particles cannot move freely along 
the stripe and the hopping in the system is still primarily 
in the vertical direction. From Fig. |l^, we can also see, 
after trying some moves, that we cannot move one par- 
ticle far away enough so as to exchange the position of 
two particles. This is a necessary condition for the boson 
and fermion spectra are identical. 

We also need to show that periodic boundary condi- 
tions do not affect the spectra, and we use again the 
idea outlined in Sec. II C. We have, as we did for the 



one-stripe- with-a- hole problem in Sec. IV, numerically 
checked a number of systems with two stripes, for exam- 
ple 4 X 10, 6 X 8, and 8 x 10, and we find that for the 
cases we checked the function a {a) is well defined. 
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FIG. 13: Two configurations of two merged stripes. Tfie 
arrows indicate the possible hops of one particle along the 
stripe. The horizontal movement is limited to the adjacent 
columns only. 





B. Stripe-stripe interaction 

^From diagonalization of two stripes with stripe length 
= 4, 6, 8, we observe that the highest-weight states in 
the ground state eigenvector have stripes far apart from 
each other, which suggests that stripes may repel. To 
study the interaction between two stripes quantitatively, 
we define the following function, 



■id) 



-^2stripes (-^a; 7 -^y ) '^-^ stripe (-^a:) 



(6.1) 



where d is the distance between adjacent stripes (here 
d — Ly/2 for two evenly spaced stripes), E2stripes the 
energy of the two-stripe system, and Estripe the energy 
of one single stripe. <j) is the energy cost per unit length 
per stripe due to stripe interaction, and it is positive 
for repelling stripes and negative for attracting stripes. 
(When we are not considering the dependence of (d) 
on Lx, we will simply use 4>{d).) In Fig. 14 we show 4>{d) 
for stripe length 4, 6, and 8. For all cases, 4>[d) is positive 
(the stripes repel) and decays as the stripe separation d 
increases. Observe that <j){d) for L^; = 8 is very much 
smaller than that for — 4, 6. The behavior of = 
8, which is not yet explained, is highly significant for 
extrapolation to the thermodynamic limit. The stripes 
with Lx — 4:,6 were so short that, from Sec. II] up to here, 



their Ly size dependence could be explained in terms of a 
particle moving in one dimension (y), completely ignoring 
their internal fluctuations. 

It is clear that as d ^ c», <j>{d) and from the 
graphs it does not decay exponentially. We try the fol- 
lowing power-law fitting function. 



AjLx) 



(6.2) 



In Fig. 15 we plot lii{4>{d)) vs ln((i) for —4,6 with 
d = 6, 8, 56, i.e., the largest lattice for i^, = 4 is 4x 112 
and for = 6 is 6 X 112. We see that the power- law 
assumption is good for — 4: with the decay exponent 
a (slope in Fig. 15) close to 2. For = 6 we observe 
that the slope approaches 2 as the stripe separation d is 
large. 
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FIG. 

(Eq. (|6.l[)) for stripe length 4, 



Stripe-stripe interaction energy cost (^(d) 
6, and 8. d is the distance 
between adjacent stripes {d = Ly/2). The lines are those 
that connect the data points. 



C. Power-law decay and stripe effective mass 



In Sec. IV B, we explained the exponential decay 



of the one-hole-with-a-stripe energy in Ly by mapping 
the stripe motion to a one-dimensional problem with 
a double-well potential. Here with two stripes in the 
system, we have shown again that stripe fluctuations 
as a function of x are limited, and we expect a one- 
dimensional potential can be sufficient in capturing the 
essential physics. 

Here we consider the two stripes as hard-core particles 
of mass TO* moving in the y direction only. This can be 
mapped to a pair of fermions on a one-dimensional ring of 
length Ly, mathematically just like the one-dimensional 
fermions of Sec. [II D| , that move in the x directicw, but 
with a completely different relation to the physics.Ea The 
wavevectors of these one-dimensional fermions are thus 
±Tr/Ly, so the total energy is 



E = 



2m* \ Ly J m*L'^, 



(6.3) 



V 



For the two-stripe problem, from curve fitting in 
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FIG. 15: log-log plot of energy cost function ^(d) (Eq. ( |6.l[ )) 
for stripe length = 4,6, with d — 6,8, ...,56. The slope 



(decay exponent a in Eq. (3.2)) is close to 2 for = 4, 
and for Lx = 6 the slope approaches 2 as stripe separation d 
increases. Linear regression is used for all L^, = 4 data, and 
for the Lx ~ 6 data, only the ten points with the largest d 
are used. 



Fig. 15, the decay exponent a is close to 2 for ~ 4. Us - 



ing Eq. ( |6.2| ) for (/)((i), the definition for in Eq. (3.1), 
and d — Ly/2 for two evenly spaced stripes, we have the 
following formula for the two-stripe interaction energy 



SL^A 



(6.4) 



where A is the factor in Eq. ( |6.2| 



Eqs. (6.3) and (S.4) give an expression for the effective 
mass of a stripe from two-stripe interaction, 



(6.5) 



The linear fitting intercept in Fi g. [Tsl gives, for ~ 4, 
A = exp(— 1.596418), and Eq. (6.5) then gives us the 
effective mass m*(4) = 1.5222. In Table ||, the 4 
stripe effective mass calculated in the one-hole-with-a- 
stripe problem is 1.3561, and in Table |l| that from the 
two-hole problem is 1.3417 for fermions and 1.3358 for 
bosons. All these results for the effective mass of a short, 
length-four stripe are comparable. 

It is clear that the shorter the stripes the better the 
one-dimensional approximating model is. For = 6, it 
can be seen in Fig. |l^ that the exponent approaches 2 in 
the large-d limit, but for d close to 56 (the largest system 
that we calculated), linear regression still gives 1.91. The 
intercept for = 6 in Fig. is not yet sufficient for us 
to use Eq. (|6.5|) to calculate the stripe effective mass. 



D. Three and four stripes 

We have also studied the interaction of three and four 
stripes with stripe length — 4. From diagonalization, 
we find that for 2, 3, and 4-stripe fermion systems we 
computed, the ground state energy always appears in the 
k = (0, 0) sector, and the highest-weight states in the 
ground state eigenvector have stripes far apart. To study 
the interaction of t hree stripes, we define, in the same 
fashion as Eq. (3.1), the energy cost per unit length per 
stripe, 



^3stripes i^-^x i -^y) ^-E^ stripe {J-'x) 

3Lx 



(6.6) 



where d — Ly/'S here. In Fig. we plot (1)3 (d) vs d and 
ln(</)3(d)) vs ln(c?). First of all, with three or more stripes, 
unlike the two-stripe case, particles have enough room to 
exchange with each other when the stripes merge, and. 
the fermion and boson energies are no longer the same.E2l 
(However, we see from the graph that the boson and 
fermion energies are only slightly different.) It is clear 
that the stripes repel and the exponent of 1^3 (d) for = 
4 is 2.004 for the boson case (and practically the same 
for the fermion case, not shown in the graph), close to 
1.983, the exponent for the two-stripe — i case (see 
Fig.|l|). 

The stripes are short enough in the x direction to be 



treated as rigid (See Sec. III). They are, in effect, three 



hardcore particles moving in one dimension. Apart from 
the fact that this one dimension runs along y, exactly 
the same problem was studied in Sec. Ill C] , where we 
recalled that it is convenient to convert such particles to 
fermions. 

Let us consider the ground state of three fermions of 
mass m* in a length-iy box. They have wavevectors 
Qy ~ and ±27r/Ly, so the total energy is 



E 




(6.7) 



For the three-stripe problem, we use the same 1/d^ rela- 
tion for (l)3{d), as in Eq. (6.2) for the two-stripe problem. 



d^ ' 



(6.8) 



where A^^L^) is a constant. Rearranging Eqs. (6.7) and 
( |6.8[ ) gives us a formula for stripe effective mass m*{Lx), 
similar to the two-stripe formula in Eq. (6.5), 



47r" 



(6.9) 



Using the linear fitting intercept in Fig. |l^, we get 
for Lj. = 4, ^3(4) = exp(- 1.394246), and then from 
Eq. (|6.9|), we get to* (4) = 1.4738, which is consistent 



with the two-stripe result 1.5222 calculated in Sec. VIC 
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ground, of the sort explained in Subsec. IIB 1. However 



e boson 

□ □ fermion 




) Boson, ln(phi)=-1 .394246-2.0041 1 9 ln(d) 



1.0 1.2 1.4 



in the smaUish systems tractable by exact diagonaliza- 
tion, the ground state would p robabl y be a stripe-array 



state, for (as we find in Sec. VII B ), the bulk energy 
density difference is quite small between the stripe-array 
state and the phase-separated state, yet the latter state 
must pay a considerable extra cost of the phase-boundary 
line tension times the droplet's perimeter. 

On the other hand, the stripe-array state in a finite 
system is extremely sensitive to boundary conditions: 
the concentration of holes may be just right to form one 
stripe across the system, but if the boundary vectors are 
both even, the system can only support an even num- 
ber of stripes. Even if the stripe-array is the correct 
thermodynamic answer at this hole filling, the best state 
available to this finite system is a separated droplet. 

We conclude that one cannot accept the behavior of 
our finite systems as a direct picture of the thermody- 
namic limit, because finite systems are dominated by 
finite-size and topological effects. The correct approach 
is to determine the equation of state for each competing 
phase by extrapolating it to the thermodynamic limit, 
and only then to perform a Maxwell construction to de- 
termine the phase stability. Note that different extrapo- 
lation schemes may be appropriate for qualitatively dif- 
ferent phases. 



FIG. 16: Stripe-stripe interaction energy cost 03 (d) 
(Eq. (6.6)) for three stripes with length 4 {Lx =4). d is 



the distance between adjacent stripes (d = Ly/3). The boson 
and fermion energies are slightly different. The log- log plot is 
shown below, showing for bosons an exponent 2.004. 



We have also calculated the energy for the smallest 
system with four stripes: the 4x12 system with 16 par- 
ticles. Here the boson energy is —10.8525 and the fermion 
energy —10.8418, both higher than the energy of four in- 
depedent stripes with length 4, —11.3137. 



VII. THE STRIPE-ARRAY 

In this paper, we have studied the limit near the half- 
filled, checkerboard state, with cases of a single stripe, 
holes-on-a-stripe, and a few stripes. As we approach in- 
termediate fillings, around 1/4, we could have two sit- 
uations: the separation of the system into particle-rich 
and hole-rich regions (sometimes argued to be quite gen- 
eral Jajp-interacting fermions with short-range interac- 
tiona^H) or else an array of stripes. How do we know 
which state is the ground state for our model? How do 
we know whether the stripe-array state is stable? 

If the true behavior, in the thermodynamic limit, is 
phase separation between the half-filled checkerboard 
state and a hole-rich liquid, then - in a sufhciently large 
finite system - the ground state would consist of one 
large (but still immobile) droplet in a checkerboard back- 



A. Stripe-array chemical potential 

We first calculate the chemical potential of the stripe- 
array, i.e., the energy per hole in creating a new stripe; 
this information will be needed for stability analysis later. 
We have prej4pusly presented this calculation in a con- 
densed form.c^ 

Say we have a L^x Ly system with p stripes stretching 
in the x direction. The number of particles per column 
is {Ly — p)/2 and the total number of particles is M = 
{Ly — p)Lx/2. Denote the total number of lattice sites 
N — LxLy, then the particle density is 



M 



P 
2L,. 



1 

2d' 



(7.1) 



where d is the separation between adjacent stripes and 
for even distributed stripe we have d = Ly/p. We then 
have d = 1/(1 - 2n). 

The energy of the stripe-array has two contributions: 
the energy from indepedent stripes and the energy due 
to stripe-stripe interaction. The energy per length (Tq 
of an infinitely long stripe is ctq = — 2/7r, which can be 
obtained from mapping the stripe to the spin- 1/2 chain.E£l 
Therefore the energy of p independent stripes is 



(7.2) 



where we have used p = Ly/d. On the other hand, the 
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energy due to stripe-stripe interaction is 



^^interaction — P0(^)-^a: 



— LxLy, 



(7.3) 



where we have used the two-stripe inter acti on energy per 
length per stripe defin ed in Eq. (3.1). (Note that 
we have seen in Sec. VI B| that (j> also depends on the 
length of the stripe L^- Fortunately, in the following, 
we will only work in the d oo limit of (j){d), which is 
independent of L^.) Combining Eq. (|7.2D and Eq. (7.3), 



we obtain the energy density of the stripe array as a 
function of particle density n, 



£sa (n) = 



indep 



interaction 



N 



= {ao + 0(rf))- = (ao + cj)idj){l - 2n). (7.4) 

Furthermore, we have seen, when we studied stripe-stripe 
interaction, (j){d) is positive and (j){d) ^ as d — > oo. So 
in the infinite-lattice limit we have 



£sain) = (Tq{1 ~ 2n). 



(7.5) 



That is to say that the chemical potential of the stripe- 
array is 



%M = _2ao = - = 1.273, 
dn TT 



(7.6) 



where the result for is used. Eq. (7.6) says that the 
slope of the stripe energy density curve is 1.273, and if 
we have an array of stripes in the system then adding a 
particle to the system will raise the energy by 1.273, or 
equivalently, adding a hole to the stripe-array will lower 
the energy by 1.273. This is an important quantity in 
the stability analysis that is to follow. 



B. Stripe-array stability 

We consider three stability conditions for the stripe- 
array agains phase separation. 



a hole to an infinitely long stripe lowers the energy by 
0.66. Here we have shown that adding a hole to a stripe- 
array lowers the energy by 1.273. The stripe-array state 
is therefore preferred. 



3. Stability condition 3: stripe- array heats phase separation 

The third stability condition comes from free energy 
analysis in statistical mechanics. Here because we are 
considering zero temperature physics and F = E — TS, 
energy is free energy. In Fig. ^ we show the two cases. 
The dashed tie-line is the Maxwell construction. It is tan- 
gent to the liquid curve (small fillings) and is connected 
to the half- filled state at n = 1/2. It represents the co- 
existence of the liquid state and the half-filled state, i.e., 
the phase separated state. According to Sec. VII A, the 
stripe-array line should have slope 1.273, and it is drawn 
from the half- filled limit. The stripe-array case is stable 
when the stripe-array line is below the dashed line, oth- 
erwise the phase separated state is stable. Therefore, to 
determine the stability of the stripe-array against phase 
separation, we need to determine the slope of the dashed 
line fi^'-^. Here we are following the notation in Ref. ^ 
and the superscript LC denotes the two states: liquid and 
CDW (charge-density- wave, i.e., half-filled state) that are 
connected by the dashed line. 
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1. Stability condition 1: stripes repel 

We have already studied stripe-stripe interaction, and 
we know that stripes repel (see Fig. ^for example). This 
is our first stability condition. 



2. Stability condition 2: two stripes heat one fat stripe 

Our second stability condition comes fro m the work on 
the one-hole- with-a-stripe problem in Sec. [VC. Here we 
want to know whether with more holes new stripes will 
form or the existing stripes just get wider (thus leading to 
phase separation). We showed in Sec. IV C| that adding 



FIG. 17: Stability of the stripe-array vs phase separation. 
Energy density E/N vs particle density n = M/N. The 
stripe-array curve has slope fj,* — 1.273. The dashed line is 
tangent to the liquid curve and is connected to the half-filled 
state. It represents the coexistence of the two phases: liquid 
and half-filled states, i.e., a phase separated state. On the 
left is the case where the stripe-array line is below the dashed 
line, and therefore the stripe-array is stable. On the right is 
the case where the dashed line is below the stripe-array line, 
i.e., the phase separated state is stable, n* is the intersection 
of the tie-line with the liquid curve, and the key quantity is 
the slope of the dashed line 

To obtain j/"'-^ , we need a fitting function for the energy 
density E/N at intermediate fillings (the liquid part). We 
use a polynomial form, up to the third order in particle 
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density n = M/N, 



8{n) 



E{n) 
N 



(7.7) 



where E{n) is the energy of the M particle system. We 
should emphasize that this fitting form Eq. (7.7) is not 
meant for the dilute limit {n —>■ oo). This is for the 
intermediate-filling, with n « 1/4. Eq. (7.7) can also be 
written as 



— ^ Ai + A2n + Asn^ 



(7.8) 



wher e E /M is the energy per particle. We will use 
Eq. (7.8) to fit diagonalization results. The slope 
can be determined by first solving for n*, the n coordi- 
nate of the intersection of the dashed tie-line with the 
liquid curve, using 



,LC 



£{n*) - _ d£{n) 



1/2 



dn 



(7.9) 



number of lattices with lattice size ranging from 25 to 
42 .E3 In Table IV we list the quadratic fitting parameters 
obtained from diagona lizat ion data. For the ten lattices 
in Fig. |l|, using Eq. ([tqI) , we obtain fj.^^ 1.254 for 



fermions and /n^*-^ = 1.304 for bosons. 



TAB LE IV: Fermion and boson fitting parameters (in 
Eq. (7.7)) and the calculated n* and ^J}"'^^ (using Eq. (7.9)), 



for all ten lattices in Fig. and for the smallest and largest 
five lattices separately. Compared with the stripe-array slope 
fi* = 1.273, the fermion stripe-array is stable against phase 
separation {fi^'^ < fi*), and the boson stripe-array is unstable 



Particle 


Lattice 


Ai A2 As 


n* m'-'^ 


SA Stable? 


Fermion 


all 10 
small 5 
large 5 


-5.454 23.597 -27.186 
-5.995 28.230 -36.825 
-5.102 20.586 -20.932 


0.251 1.254 
0.255 1.240 
0.250 1.264 


Stable 


Boson 


all 10 
small 5 
large 5 


-6.138 27.839 -33.950 
-6.378 29.929 -38.420 
-5.994 26.580 -31.251 


0.233 1.304 
0.231 1.301 
0.233 1.307 


Unstable 





E/M=-6.138+27.839n-33.950n'n 



-2.6 I 

0.15 



FIG. 18: Fermion (top) and boson (bottom) E/M vs n = 
M/N fit for intermediate fillings (0.19 < n < 0.31). The 
fermion data have larger spread than the boson data because 
of the fermion shell effect 

In Fig. |lj, we plot, for fermions and bosons, E/M vs 
n at the intermediate fillings (0.19 < n < 0.31) for a 



In Table IV, we have also included calculations for 
the five smallest and five largest lattices in Fig. |lj sep- 
arately. We find that f/"'^ for both boson and fermion 
increases with increasing lattice size. For the boson case, 
this means that for large lattices, the stripe-array is not 
stable against phase separation (/i'^^ > fi* = 1.273). For 
the fermion case, the result is stable (/i'"*"^ < /i*) for the 
lattices that we have studied, but n^'-^ is very close to the 
stability limit 1.273. The fermion stripe- array is possibly 
stable. 

In our earlier publication, Ref . we used the same 
polynomial fitting function, Eq. (7.7), but we fixed Ai = 



—4, corresponding to the energy of an noninteracting par- 
ticle, and our systems there went from 20 sites (4 x 5) 
to 36 sites (6 x 6). The analysis in this section and 
the results in Table |^ are obtained using the three- 
parameter (^1,^12,^3) fit for bigger lattices, from 25 
sites (5 X 5) to 42 sites (6 x 7). The conclusion in Ref. |35| 
was fi^'^ = 1.25(2) for fermions, and fi^'^ = 1.33(2) for 
bosons. 

R ef. p4 simulated the boson case of our model 
Eq. ( |1.1D , using Quantum Monte Carlo, but with < c». 
Their phase diagram, Fig. 11 of Ref. p4l shows a first- 
order phase transition between the (tt, tt) solid [i.e. the 
checkerboard phase] and a superfluid [the hole-rich phase 
with n < 0.25. Our results correspond to a single point 
{t/V,fi/V) = (0,0) on their plot; however, tijeir phase 
boundary approaches that point with a slopeE3 2.28, i.e. 
at ^/t — 2.28 inXhe V 00 limit, or fi/t — 2.4 in a re- 
lated simulationE3. This is much larger than our result 
fJ-LG ~ 1.33. That is puzzling, since it would imply either 
that the liquid state (with density n « 0.25) has much 
lower energy than we found, or else that a state exists 
at n « 0.4 (in these finite systems) having much lower 
energy than the stripe states we studied. 

The evidence in Ref. |63 for phase separation rests 
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on histograms such as their Fig. 10 which exhibits (for 
V = 2.86 in an 8 X 8 system) a bimodal probabihty dis- 
tribution for the density n, in a grand canonical ensem- 
ble, with peaks aX n — 0.5 (the checkerboard state) and 
n « 0.39, the latter being interpreted as a hquid state. 
For our V = oo model, a similar plot at low temperature 
would show peaks at n = 0.5 and at n = 0.375. The latter 
peak represents a two-stripe state with 8 holes. Its energy 



-8a/1 + 1/V2 w -10.45 (using Eq. ( |3^ ), and ne- 



IS 

glecting the stripe-stripe interaction in light of Fig. 14). 
Our point is that examination of the most probable con- 
figurations is needed in order to judge whether the co- 
existing state contains stripes or liquid droplets. Indeed, 
the discrete coexistence between n = 0.5 and n = 0.375 
is an artifact of the periodic boundary conditions - in 
the thermodynamic limit, within the stripe-array phase, 
density n varies continuously with the chemical potential. 



VIII. CONCLUSION 



A. Summary 

In this paper we studied systematically a two dimen- 
sional model of strongly-interacting spinless fermions and 
hardcore bosons on the square lattice Eq. (1.1) near the 
half-filled limit. We considered an extended, quantum- 
fluctuating object that is natural in our model at the 
dense limit-the stripe. After introducing our model 
(Sec. I) and our diagonalization program (Sec. ||), we 
made a detailed study of the problems of a single stripe 
(Sec. 01), up to two holes on a stripe (Sees. |^ and 0), 



stripe-stripe interaction of two, three, and four stripes 
(Sec. |vi|) , and finally the stability of an array of stripes 



(Sec. Vn). Our theoretical studies were aided with sub- 
stantial amount of results obtained from exact diagonal- 
ization. 

As is well-known, quantum systems in two dimensions 
are much more difficult to study than those in one di- 
mension. Because of the presence of strong repulsion, at 
the dense limit of our model, motion is severely limited. 
We showed that for the problems of a single stripe, holes 
on a stripe and a few stripes, analytical results can be 
obtained by mapping the two-dimensional problems to 
one-dimensioinal ones, which can in turn be studied us- 
ing the analytical tools available in one dimension. We 
showed in great detail the intricacies of each map, includ- 
ing for example, the many-to-one problem, the resulting 
phase factors, and the role of boson and fermion statis- 
tics. 

Another theme in this paper is the comparison between 
the boson and fermion problems. At the center of this 
comparison is the state graph introduced in Sec. || which 
describes graphically the relationship among the basis 
states. For the problems of a single stripe, holes on a 
stripe, and a few stripes, we investigated the effect of par- 
ticle statistics on eigenenergy, with the help of the state 



graph and our computer program. And for the problem 
of stripe-array stability, we obtain from diagonalization 
results an interesting conclusion that the boson stripe- 
array is not stable and the fermion stripe-array is very 
close to the stability limit and is possibly stable. 

This paper taken as a whole is an attempt to under- 
stand stripes as macroscopic objects arising from micro- 
scopic components (the spinless fermion and hardcore bo- 
son particles in our model) . The underlying physics is the 
quantum mechanics of the component particles that is de- 
scribed by a many-particle Schrodinger equation (which 
we solve by exact diagonalization). Stripes can perhaps 
be called an emergent phenomenon and are results of col- 
lective motions of many microscopic particles (the origi- 
nal particles or holes). They can perhaps be considered 
in the same fashion as phonons arising from collective 
lattice vibrations. And like phonons, the stripes in our 
problem have a life of their own, and we have tried, in 
this paper, to understand the new physics they bring. 

This paper is one of the first using exact diagonal- 
ization to study stripes. This is possible because of the 
signi fica nt Hilbert space reduction afforded by our model 
Eq. ( p_.l| ) , spinless fermions and hardcore bosons with in- 
finite nearest-neighbor repulsion. We believe that this is 
one of the simplest microscopic lattice model with which 
to study stripes from the underlying particle dynamics. 
And as we commented in our other paper on the dilute 
limit of this model Jj this model deserves to be better 
known and better studied. 



B. Future directions 

There are a number of directions in which the work of 
this paper might be extended. For example, our code 
handles quite arbitrary periodic boundary conditions, 
not necessarily rectangular or square. Our key use of this 
feature was to accurately measure the equation of state 
for the strongl y inter acting liquid at density n « 0.25, 
as used in Sec. VII B to decide the stripe-array stability. 
By including a wide variety of differently shaped fermion 
systems in Fig. |l^ (top plot), in effect we average out 
the "noise" (which is due to fermion shell effectsu). An 
alternative way to do the same thing, which we have not 
implemented in the present work, would be to impose a 
phase factor acpass the boundary conditions and average 
over all phases ED 

A boundary condition could be used to investigat e the 
90° bend of a stripe, which was discussed in Sec. VC. 
These would force a diagonal stripe, as in Fig. |ll|, but 
additional holes would be added to the stripe which are 
expected to condense and form a segment at right angles. 
That geometry would permit estimation of the energy 
cost of the 90° bend. 

This study has focused entirely on eigenenergies; since 
our diagonalizations provide the wavefunctions too, we 
could have computed a variety of informative expecta- 
tions, such as correlation functions. We did make quali- 
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tative observations using snapshots of the high- weighted 
configurations. In our diagonalizations with two or more 
stripes (Sec. VI D), which started from a state with 
stripes all merged, the highest-weight state in the ground 
state eigenvector has stripes far apart from each other, 
and nearly equally spaced, as expected from the one- 
dimensional approximate wavefunctions of Sec. Cal- 
culation of the probability distribution for the stripe sep- 
aration would provide a quantitative test of those wave- 
functions. Similar calculations would give a more direct 
check of the exponential decay of the hole probability as 
a function of its distance from a stripe, as predicted in 
Sec. |VB| ; and would immediately reveal the attra ctive or 
repulsive tendency of two holes on a stripe (Sec. VC). 

We studied the case of one stripe with two holes; the 
case of two stripes with one hole is also worth investigat- 
ing. The finite-size dependence of the extra energy due 
to the hole would shed light on the processes by which 
the hole is transferred between stripes when they touch. 
Such processes are critical to the transport properties of 
the stripe array, as considered in Refs. M in the spinfull 
case. 

Up to here, we discussed further kinds of measurements 
which could be made on the same model system; varia- 
tions are also possible of the model itself. The most ob- 
vious of these is to treat V < oo. A large but finite V 
may be handled in the spirit of the t-J: the same highly 
restricted basis states (which facilitated our diagonaliza- 
tions) are retained, but new hopping terms along [2,0] 
and [1,1] type vectors appear, with amplitudes — t^/F, 
wherever the intervening site would be forbidden by the 
nearest-neighbor exclusion. It is not obvious whether this 
tends to stabilize or to destabilize the stripe array. The 
inclusion of another form of correlated hopping (Ref. |3^, 
Fig. 3(c)) can certainly stabilize or destabilize the stripe 
array depending on its sign.ta 

Finally, the same code is adaptable, almost without 
modifications, to the triangular lattice. That may model 
•^He or iJJe-atoms adsorbed on graphite or on carbon nan- 
otubesJlj'Ej which implement a periodic boundary con- 
dition in one direction. 
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APPENDIX A: COUNTING BASIS STATES FOR 
THE PROGRAM 

We find the number of basis states Af is maximum 
for fillings n = M/N !^ 0.25, close to the filhng n* at 



which, in the thermodynamic limit, (see Sec. VII B) a 
hole-rich liquid coexists with the half-filled state or the 
(nearly-half-filled) stripe array state. This count may be 
estimated by adapting Pauling's trick for the ice model 
entropyES. The number of ways merely to distribute M 
particles over N sites, uncontrained, is (^) implying the 
usual entropy —n Inn — (1 — n) ln(l — n) as TV — > oo; this 
must be corrected to account for the constraint of no 
nearest neighbors. One chosen particle has a probability 
(1 — n)** to be free of nearest neighbors. If we pretend 
this event is independent as each particle on (say) the 
even lattice is chosen in turn, then the M/2 power of 
this probability is the chance for the whole configuration 
to be valid. The net entropy is estimated as — nlnn — 
(1 — 3n) ln(l — n); this attains its maximum 0.42 around 
n — 0.24. implying the leading dependence Af ^ 1.54^ 
for the U = oo spinless case, compared to Af ^ 4^ in 
the Hubbard model case. In fact, for a 7 x 7 lattice, the 
block with M = 11 particles has the largest (translation- 
reduced) matrix dimension, Af = 1906532. 

Let Afjl{Lx, Ly) be the number of basis states (per k 
vector) in a configuration with s stripes and h holes on 
it. This grows exponentially with L^; but comparatively 
slowly with Ly, so we can handle quite large Ly, and 
moderately large L^. 

In the c ase of one stripe, the map to an XX spin chain 
(Sec. Ill Bp makes clear that, independent of Lj,, the basis 



has ~ L^'\lJ2) states (where (^) = A\/[B\{A^ By.]). 

Asymptotically Af^ (L^, Ly) - 2^- /LJ^^. 

The number of states with two stripes may be esti- 
mated by placing them independently - the correction 
from disallowing overlaps is subdominant. That is the 
square of the one-stripe count, multiplied by L^Ly for 
the possible vectors offsetting one stripe relative to the 



other, 



AfiiL, 



Ly) 



3/2 



A^-Ly/Ll 



For one stripe with one hole, the basis states were 
counted for Ly < 25 and < 10. Empirically, 



AAi(L„L,) = - 



1/ 



2 \L^/2 



{Ly-l) + B{L,] 



(Al) 



exactly, where B(4) = -2, B(6) = -4, B(8) = -8, 
S(10) — —14; for Ly — 7, this appears to follow with 
S(12) = -12, B{U) = 66, B(16) = 572. The linear in- 
crease with Ly is due to the states with the hole away 
from the stripe. 

For a stripe with two holes, Eq. ( Al) gets multiplied by 
an additional factor L^Ly, for the possible offset of the 
second hole relative to the first. Empirically, A/'j^ (6, Ly) — 
IbLy^ — 4:7Ly + 4 (valid for Ly > L^); this is inferred 
from Ly < 19. For Ly — 7, we found roughly Af{Lx, 7) ~ 



3.5L^^/^2'^-, for L^ up to 13. 
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